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Abstract 

Given k pairs of terminals {(si,fi),..., (s/t, ip} in a graph G, the min-sum k vertex-disjoint 
paths problem is to find a collection {Qi, Q2, ■ ■ ■, Qk} of vertex-disjoint paths with minimum 
total length, where Qi is an Sj-to-fi path between .s, and ti. We consider the problem in 
planar graphs, where little is known about computational tractability, even in restricted cases. 
Kobayashi and Sommer propose a polynomial-time algorithm for k < 3 in undirected planar 
graphs assuming all terminals are adjacent to at most two faces. Colin de Verdiere and Schrijver 
give a polynomial-time algorithm when all the sources are on the boundary of one face and all 
the sinks are on the boundary of another face and ask about the existence of a polynomial-time 
algorithm provided all terminals are on a common face. 

We make progress toward Colin de Verdiere and Schrijver’s open question by giving an 
0(kn 5 ) time algorithm for undirected planar graphs when {(si, U),..., (sfc,ip} are in counter¬ 
clockwise order on a common face. 
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1 Introduction 


Given k pairs of terminals {(si,ii), ..., (s k ,t k )}, the k vertex-disjoint paths problem asks for a set 
of k disjoint paths {Qi,Q 2 , ■ ■ ■ ,Qk}, in which Qi is a path between s, and t, for all 1 < i < k. 
As a special case of the multi-commodity flow problem, computing vertex disjoint paths has found 
several applications, for example in VLSI design |KvL84j, or network routing [ORS93, SM05|. It is 
one of Karp’s NP-hard problems [Kar74 even for undirected planar graphs if k is part of the input 
[MP93 . However, there are polynomial time algorithms if A: is a constant for general undirected 

In general directed graphs, the A:-vertex-disjoint paths problem is NP-hard 
but is fixed parameter tractable with respect to parameter k in directed 


graphs RS95 KW10 


even for k = 2 


FHW80 


planar graphs [Sch94 , CMPP13 


Surprisingly, much less is known for the optimization variant of the problem, minimum-sum 
k vertex-disjoint paths problem (k-min-sum), where a set of disjoint paths with minimum total 
length is desired. For example, the 2-min-sum problem and the 4-min-sum problem are open 
in directed and undirected planar graphs, respectively, even when the terminals are on a common 
face; neither polynomial-time algorithms nor hardness results are known for these problems |KS10 


Bjorklund and Husfeldt gave a randomized polynomial time algorithm for the min-sum two vertex- 
disjoint paths problem in general undirected graphs BH14|. Kobayashi and Sommer provide a 
comprehensive list of similar open problems (Table 2 | KS10]). 


One of a few results in this context is due to Colin de Verdiere and Schrijver |VS111: a polyno¬ 


mial time algorithm for the A-rniri-surri problem in a (directed or undirected) planar graph, given 
all sources are on one face and all sinks are on another face [V SlT] , In the same paper, they 
ask about the existence of a polynomial time algorithm provided all the terminals (sources and 
sinks) are on a common face. If the sources and sinks are ordered so that they are in the order 
si, S 2 , ■ ■ ■, s k ,t k ,t k -i, ■ ■ ■ ,t i around the boundary, then the A;-min-sum problem can be solved by 
finding a min-cost flow from si, S 2 , ■ ■ ■, s k to t k ,tk-i, ■ ■ ■, t\. For A < 3 in undirected planar graphs 
with the terminals in arbitrary order around the common face, Kobayashi and Sommer give an 

KSloF 


0(n 4 log n) algorithm 


In this paper, we give the first polynomial-time algorithm for an 


arbitrary number of terminals on the boundary of a common face, which we call F, so long as the 
terminals alternate along the boundary. Formally, we prove: 

Theorem 1.1. There exists an 0{kn 5 ) time algorithm to solve the k-min-sum problem, provided 
that the terminals si, t\, S 2 , t 2 , ■ ■ ■, s k , t k are in counter-clockwise order on the boundary of the graph. 


Definitions and assumptions We use standard notation for graphs and planar graphs. For 
simplicity, we assume that the terminal vertices are distinct. One could imagine allowing t* = Sj+i; 
our algorithm can be easily modified to handle this case. We also assume that the shortest path 
between any two vertices of the input graph is unique as it significantly simplifies the presentation 


of our result; this assumption can be enforced using a perturbation technique MVV87 . 


2 Preliminaries 

Walks and paths. Let G = (V,E) be a graph, and let LA be a subgraph of G. We use V{H) 
and E(H) to denote the vertex set and the edge set of AA, respectively. For U C V, we use G[U] 
to denote the subgraph of G induced by U, whose vertex set is U and whose edge set is all edges 
of G with both endpoints in U. A walk W = (v\,V 2 , ■ ■ ■ ,v k ) in G is a sequence of vertices such 

1 Kobayashi and Sommer also describe algorithms for the case where terminals are on two different faces, and 
k = 3. 
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that for all 1 < i < h, we have {vi,Vi + \) G E. For any Vi,Vj G W, W[vi,vj\ is the (sub-)walk 
(vi,Vi+ 1 ,... ,Vj). A path is a walk with no repeated vertices. Sometimes, we view a walk as its set 
of edges, and use set operations on walks. For example, given two walks W\ and W 2 , we define 
W = W\ 0 W 2 to be their symmetric difference when it is clear from the context that W is a walk. 
Given a length function £ : E —> V, the length of W is denoted by £{W), and it is ^• ! = T 1 1 &{ v ii v i+i)- 

Planarity. An embedding of a graph G = (V, E) into the Euclidean plane is a mapping of vertices 
of G into different points of M 2 , and edges of G into internally disjoint simple curves in M 2 such 
that the endpoints of the image of ( u, v) G E are the images of vertices u G V and vGh If such 
an embedding exists then G is a planar graph, A plane graph is a planar graph and an embedding 
of it. The faces of a plane graph G are the maximal connected components of M 2 that are disjoint 
from the image of G. If G is connected it contains only one unbounded face. This is called the 
outer face of G and we denote the boundary of G by dG. Let W be a walk and let W° be the set 
of edges that are used an odd number of times in W. W° is a collection of simple cycles. For any 
point p G R 2 , we say that p is inside W if p is on the image of W in the plane, or c is inside an odd 
number of cycles of W°. When it is clear from the context, we use the same notation to refer to 
the subgraph composed of the vertices and edges whose images are completely inside W. 

3 Structural properties 

In this section, we present fundamental properties of the optimum solution that we exploit in our 
algorithm. To simplify the exposition, we search for pairwise disjoint walks rather than simple 
paths and refer to a set of pairwise disjoint walks connecting corresponding pair of terminals as a 
feasible solution. Indeed, in an optimal solution, the walks are simple paths. 

Let {Q 1 , Q 2 , ..., Qk} be an optimal solution, where Qi is a Sj-to-L path and let {Pi, P 2 ,..., P*,} 
be the set of shortest paths, where P* is the si-to-ti shortest path. These shortest paths together 
with the boundary of the graph, dG, define internally disjoint regions of the plane. Specifically, 
we define Ri to be the subset of the plane bounded by the cycle P* U C t , where C, is the Si-to-ti 
subpath of 8G that does not contain other terminal vertices. The following lemmas constrain the 
behavior of the optimal paths. 

Lemma 3.1. For all 1 < i < k, the path Qi is inside R, % . 

Proof. Suppose, to derive a contradiction, that Qi <£. Ri. So, there is a vertex v G Qf\Ri. Let p v 
be the maximal path containing v that is internally disjoint from Pj. Let ( x , y ) be endpoints of p v , 
and observe that x,y G Pi- Also, by uniqueness of shortest paths, we have £{P t \jx, y\) < £(Qi[x,y]). 
Thus, Ql = Qi® Qi[x, y] 0 P t [x, y] is shorter than Qi. 

It remains to show that Q' t n Qj = 0 for all j 7 ^ i, 1 < j < k. But, by the construction, Q\ 
is inside Ci U Qi, and all terminals other than si and U are outside C t U Qi. So, by Jordan curve 
theorem, for any Qj to intersect Q\, it has to intersect Qi, too. Thus, {Qi,..., Qk}\Qi U Q\ is a 
shorter solution than the optimum, which is a contradiction. □ 

We take the vertices of P t and Qi to be ordered along these paths from Sj to Q. 

Lemma 3.2. For u,»G Qi H P,: ; u precedes v in Pi if and only if u precedes v in Qi. 

Proof. Suppose, to derive a contradiction, that u and v have different orders on Pi[si,U] and 
Qi[si,ti], and assume, without loss of generality, that u precedes v in P,[s*, tf , but v precedes u in 
Qi[si,ti\. So, Qi[si,U] can be decomposed into three subpaths (1) 71 is a Si-to-v path, (2) 72 is a 
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v-to-u path, and (3) 73 is a it-to-f,; path. If 71 contains u then Qi is not a simple path, visiting u at 
least twice. Otherwise, The Jordan curve theorem implies that 72 or 73 must intersect 71 . Again 
Qi is not simple, so, it is not a walk in the optimum solution. □ 



Figure 1: (left) A 4-min sum instance; regions are shaded and borders are green, (right) A feasible 
solution; Type I and Type II subpaths are blue and red, respectively. 


We call Ri n Rj the border of R{ and Rj and denote it Bij. Note that a border can be a 
single vertex. Since we assume shortest paths are unique, Bij is a single (shortest) path. Figure [I] 
illustrates borders for a 4-min-sum instance. The following lemma bounds the total number of 
borders. 


Lemma 3.3. There are 0(k ) border paths. 

Proof. Let 1Z = be the graph whose vertices correspond to regions in G, and there 

is an edge between two vertices of Vn if their corresponding regions in G share a border. Let 
H = G[dF U Pi U ... U Pfc], and observe that IZ is a subgraph of the planar dual of H. Thus, IZ is 
planar. 

Since there is a bijection between Vjz and the set of regions of G, we have \ Vn\ = k. Additionally, 
there is a bijection between E-ji and borders in G. Because, 1Z is planar we conclude that the number 
of borders in G is \En\ = 0(Vjz) = 0(k). □ 


Consider a region Ri and consider the borders along P t , B % i x , P v ,; 2 ,..., Observe that the 
intersections of the regions R . lx , Pj 2 ,..., Ri t with 8G must be in a clockwise order. Let n,..., n 
be the subsequence of i\,... of indices to regions that intersect Qi. For j e {ti,..., ig}, let Xj 
and yj be the first and last vertex of Qi in Bij. Additionally, define yo = Si and 27+1 = ti. We 
partition Qi into a collection of subpaths of two types as follows. 


Type I : For h = 0,... ,£, Qi[yt, Xh+ 1 ] is a Type I subpath in region Ri. 

Type II : For h = 1,... — 1, Qi[xh, yh ] is a Type II subpath in region Ri. We say that Qi[xh, yh} 
is on the border Bij containing Xh and yh- 


By this definition, all Type I paths are internally disjoint from all borders. By Lemma 3.2 each 
Type II path is internally disjoint from all borders except possibly the border that contains its 
endpoints, with which it may have several intersecting points. See Figure [l] for an illustration of 
Type I and II paths. 

The following lemma demonstrates a key property of Type I paths, implying that (given their 
endpoints) they can be computed efficiently via a shortest path computation: 


Lemma 3.4. Let a be a Type I subpath in region Ri. Then a is the shortest path between its 
endpoints in Ri that is internally disjoint from all borders. 
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Proof. Let u, v E V (G) be the endpoints of a , and let a' be the shortest u-to-v path in Ri that is 
internally disjoint from all borders. Also, let Q\ = Qi®a®a'. The path a' has the same endpoints 
as a, and it is internally disjoint from all Rj if j 7 ^ i. Also, by Lemma 3.1, for each 1 < j < k, Qj is 
inside Rj. Therefore, a! is disjoint from Qj for all 1 < j < k and j 7 ^ i. Thus, {Q 1 ,..., Qk}\Qi U Q[ 
is a set of k pairwise disjoint walks, with total length OPT — (.{a) + £{a!) where OPT is the total 
length of the optimal paths. So, OPT — £{a) + £{a!) > OPT , which implies £{a') > £(a). Therefore 
a must be a shortest (u,v) path in Ri that avoids boundaries. In fact, uniqueness of shortest paths 
implies a = a'. □ 


A Type II path has a similar property if it is the only Type II path on the border that contains 
its endpoints. The proof of the following lemma is almost exactly the same as the previous proof. 

Lemma 3.5. Let f3 be a Type II subpath in region Ri on border Bij. Suppose there is no Type II 
path on Bij inside Rj. Then [3 is the subpath of Bij between its endpoints. 

The following lemma reveals a relatively more sophisticated structural property of Type II paths 
on shared borders. 


Lemma 3.6. Let (3 and 7 be Type II subpaths in Ri and Rj on Bij, respectively, let Xi and yi be 
the endpoints of f3, and let Xj and yj be the endpoints of 7 . Then, {/9, 7 } is the pair of paths with 
minimum total length with the following properties: 

(1) (3 is an Xi-to-yi path inside Ri, and it is internally disjoint from all borders except possibly 
33ij [•£*) Vi] ■ 

(2) 7 is an Xj-to-yj path inside Rj, and it is internally disjoint from all borders except possibly 

Bij[xj,yj\. 

Proof. Properties (1) and (2) of (3 and 7 are implied by the definition of Type II paths and because 
they are internally disjoint from all borders except Bij. It remains to show that their total length 
is minimum. 

Let (/ 3',y) be the pair of paths with minimum total length that has properties (1) and (2). Let 
Q'i = Qi® (3 ® /3', and let Q' rj = Qj® 7 © 7 '. 

By construction, f3' is internally disjoint from all borders except (possibly) Bij[xt, yi\. Addition¬ 
ally, the endpoints of j3' are the same as the endpoints of (3. Therefore, intersection points of (3' with 
borders of Ri that are not in f3 are all in Bij[xi, yf\. Similarly, intersection points of 7 ' with borders 
of Rj that are not in j3' are all on Bij\xj,yj\. It immediately follows that Q\ and Qh are disjoint 
for any ft £ {1,2,..., fe}\{j}. Similarly, Q' rj and Qh are disjoint for any h 6 {1, 2,..., fc}\{i}. 

Furthermore, Q\ and Q'- are disjoint by their construction. Thus, {Qi, ■ ■ ■, Qk}\{Qi, Qj} U 
{Qii Qj} is a set of k pairwise disjoint walks connecting the terminals. The total length of this new 
set is OPT — £(/3) — £( 7 ) + t{f3') + 1{ 7 ') > OPT. Consequently, £((3) + £( 7 ) < £(/3') + £( 7 '), in fact, 
£Q3) + £( 7 ) = £(/3') + £(y)- Thus, (a, (3) are a pair of path with minimum total length that has 
properties ( 1 ) and (2). □ 


4 Algorithmic toolbox 

In this section, we describe algorithms to compute paths of Type I and II for given endpoints. 
These algorithms are key ingredients of our strongly polynomial time algorithm described in the 
next section. More directly, they imply an n°^ time algorithm via enumerating the endpoints, 
which is sketched at the end of this section. 
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Each Type I path can be computed in linear time using the algorithms of Henzinger et al. [HKRS97 


they can also be computed in bulk in 0 (n log n) time using the multiple-source shortest path al¬ 
gorithm of Klein |Kle05| (although other parts of our algorithms dominate the shortest path com¬ 
putation). Similarly, a Type II path on a border Bij can be computed in linear time provided it 
is the only path on B tJ . Computing pairs of Type II paths on a shared border is slightly more 
challenging. To achieve this, our algorithm reduces this problem into a 2-min sum problem that can 
be solved in linear time via a reduction to the minimum-cost flow problem. The following lemma 
is implicit in the paper of Kobayashi and Sommer KSlOj. 


Lemma 4.1. There exists a linear time algorithm to solve the 2-min sum problem on an undirected 
planar graph, provided the terminals are on the outer face. 


Proof. Consider an instance of the 2-min sum problem with terminals {(si, t\), (s 2 , £ 2 )}- This 
problem reduces to a minimum-cost flow problem as follows. Because of the symmetry for terminals 
in undirected graphs, we can assume that si and S 2 are next to each other on the outer face: there 
is a si-to-S 2 path on the boundary of the outer face that does not intersect 1 , ^ 2 }- We make 
the graph directed by replacing each undirected edge {u, u} with edges (u,v) and (v,u). For edge 
(u,v) in the directed graph, we assign its length using length function t(u,v). For every vertex v 
in G, we split it into two vertices v\ and V 2 and connect them with a zero length edge that has 
unit capacity. For each edge (u,v), we connect u to ui, and for each edge (v,u), we connect V 2 to 

u. We introduce a dummy source vertex d\, with two edges (di,si), and (di,S 2 ) of unit capacity 
and zero length. Also, we introduce a dummy sink vertex d 2 , with edges (t 1 , c^ 2 ) and (t 2 , ^ 2 ) of unit 
capacity and zero length. We assign capacity one to all other edges. The lengths of the other edges 
are specified by the length function i of G. Since, si and s 2 (also t\ and £ 2 ) are next to each other, 
the graph remains planar after adding u, v, and their incident edges. Now, it is straight forward 
to see that our 2 -min sum problem is equivalent to a minimum cost u-to-v flow problem of value 
2. This minimum cost flow problem in turn reduces to two shortest path computations from u to 

v, which can be done in linear time (HKRS97]. □ 


We reduce the computation of Type II paths to 2-min sum. The following lemma is a slightly 
stronger form of this reduction, which finds application in our strongly polynomial time algorithm. 


Lemma 4.2. Let Ri and Rj be two regions with border B 7 ] and let Xi,yi E P t and Xj,yj E Pj. A 
pair of paths (/3, 7 ) with total minimum length and with the following properties can be computed 
in linear time. 


1. (5 is an Xi-to-yi path inside Ri, and it is internally disjoint from all borders except possibly 

P{ [x.i, yi] n Bi j. 

2. 7 is an Xj-to-yj path inside Rj, and it is internally disjoint from all borders except possibly 
Pj [ x j 5 yj] n Bij. 

Proof. Let the graph H be the induced subgraph by the vertices of G inside Ri U Rj. We obtain 
H' from H by performing the following operations: 

1. deleting all vertices of H that belong to borders other than Bij, 

2. deleting all edges in Ri that are incident to Bij but not incident to Vi] H Bi,ji and 

3. deleting all edges in Rj that are incident to Bij but not incident to Pj[xj,yj\ n Bij. 
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Figure 2: (left) The graph H induced by vertices in Ri U Rj] regions are shaded and borders are 


green, (right) The graph H' defined in the proof of Lemma 4.2 


Note that H' is not necessarily connected. For an illustration of H and H' , see Figure [2j 

Since j3 and 7 are intact in H' , they can be computed in H' instead of G. Furthermore, observe 
that Xi, yi,Xj,yj are on the boundary of H' . If H' is disconnected, then /3 and 7 are a shortest paths 
in their own connected components. So they can be computed in linear time using the algorithm 
of Henzinger et al. [HKRS97 


So, suppose H' is connected, and observe that Xi,yi,Xj,yj are on the boundary of the outer 
face of H'. By Lemma 4.1, there is a linear time algorithm to compute a pair of disjoint paths of 
minimum total length between corresponding terminals. 

Let j3' and 7 ' be Xj-to-y* path and Xj-to-yj path computed via solving a 2-min sum problem. 
It remains to prove that f3' and 7 ' have properties (1) and (2) of the lemma. That is /3 7 € H' n Ri, 
and y £ H' n Rj. This can be done through a replacing path argument similar to Lemma 3.1 
Bij (so, any subpath of it) is a shortest path. 


as 

□. 


4.1 An n°^ time algorithm 

The properties of Type I and II paths imply a naive n° ^ time algorithm, which we sketch here. 
An optimal solution defines the endpoints of Type I and Type II paths, so we can simply enumerate 
over which borders contain endpoints of Type I and II paths and then enumerate over the choices 
of the endpoints. Consequently, there are zero, two, or four (not necessarily distinct) endpoints of 
Type I and II paths on Bi or 

'mj)\ + fm,i) 


1 + 


possibilities, which is 0(n 4 ) since £(Bij) = 0(n). Since there are 0(k ) borders (Lemma 3.3), there 
are n°^ endpoints to guess. Given the set of endpoints, we compute a feasible solution composed 
of the described Type I and II paths or determines that no such solution exists. Since Type I and II 
paths can be computed in polynomial time, the overall algorithm runs in n°^ time. 


5 A fully polynomial time algorithm 

We give an 0(kn 5 )-time algorithm via dynamic programming over the regions. For two regions Ri 
and Rj that have a shared border Bij, Ri and Rj separate the terminal pairs into two sets: those 
terminals S£, te for l = i + 1,. .. ,j — 1 and s m , t m for mn = j + 1,..., k, 1,. .., i — 1 (for i < j). Any 
S£-to-t£ path that is in region Ri cannot touch any s m -to-i m path that is in region R m since Re and 
R m are vertex disjoint. Therefore any influence the S£-to-te path has on the s m -to-t m path occurs 
indirectly through the Si-to-ti and Sj-to-tj paths. Our dynamic program is indexed by the shared 
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borders B,j and pairs of vertices on (a subpath of) P t and (a subpath of) Pj\ the vertices on Pi 
and Pj will indicate a last point on the boundary of Ri and Rj that a (partial) feasible solution 
uses. 

We use a tree to order the shared borders for processing by the dynamic program. Since there 
are 0{k) borders (Lemma 3.3), the dynamic programming table has 0(kn 2 ) entries. We formally 
define the dynamic programming table below and show how to compute each entry in 0(n 3 ) time. 


5.1 Dynamic programming tree 

Let 1Z = {Ri}i = i and B be the set of all borders between all pairs of regions. We assume, without 
loss of generality, that 1Z is connected, otherwise we split the problem into independent subproblems, 
one for each connected component of 7Z. 

We define a graph T (that we will argue is a tree) whose edges are the shared borders B between 
the regions 1Z. Two distinct borders Bjj and B^n are incident in this graph if there is an endpoint 
x of Bj -j and y of B/ that are connected by an x-to-y curve in M 2 \ F that does not touch any 
region 7 Z except at its endpoints x and y, see Figure [3j Note that this curve may be trivial (i.e. 
x = y). The vertices of T (in a non-degenerate instance) correspond one-to-one with components 
of M 2 \ (F U 1Z) (plus some additional trivial components if three or more regions intersect at a 
point), or non-regions. The edges of T cannot form a cycle, since by the Jordan Curve Theorem 
this would define a disk that is disjoint from dF ; an edge Bjj in the cycle bounds two regions, one 
of which would be contained by the disk, contradicting that each region shares a boundary with 
dF. Therefore T is indeed a tree. We use an embedding of T that is derived naturally from the 
embedding of G according to this construction. We use this tree to guide the dynamic program. 




Figure 3: (left) Thick green segments are borders, thin green curves are in M 2 \ (FUlZ) connecting 
borders that are incident in T. (right) The directed tree T used for dynamic programming. 


By the correspondence of the vertices of T to non-regions, we have: 

Observation 5.1. The borders Bj j. Bn,... along a given region Ri form a path in T. The order 
of the borders from Si to ti along Pi is the same as in the path in T. 


Consider two edges Bij and B^j that are incident to the same vertex v of T and that are 
consecutive in a cyclic order of the edges incident to v in X”s embedding. By Observation 5.1 and 
the embedding of T, there is a labeling of i,j , h,£ such that: 


Observation 5.2. Either j = h or Bij is the last border of Rj along Pj and Bh^ is the first border 
of Rh along P h . 


Root T at an arbitrary leaf. Since 1Z is connected, the leaf of T is non-trivial; that is, it has an 
edge Bij incident to it. By Observation 5.1, Bij is (w.l.o.g.) the last border of Ri along Pi and 
the first border of Rj along Pj. By the correspondence of the vertices of T to non-regions, and the 
connectivity of 7 Z, either j = * + lori = l and j = k. For ease of notation, we assume that the 
terminals are numbered so that i = 1 and j = k. We get: 





Observation 5.3. Every non-root leaf edge of T corresponds to a border Bjj + \. 

We consider the borders to be both paths in G and edges in T. In T we orient the borders 
toward the root. In G, this gives a well defined start aij and endpoint bij of the corresponding 
path Bij (note that ajj = bij is possible). By our choice of terminal numbering and orientation of 
the edges of T, from Sj to L along P t , b t j is visited before ajj, and from Sj to tj along Pj , ajj is 
visited before bij. 


5.2 Dynamic programming table 

We populate a dynamic programming table Djj for each border Bij. Dij is indexed by two vertices 
x and y: a: is a vertex of Pj,[bij. tf and y is a vertex of Pj[sj,bij\. D l , ] [x. y\ is defined to be the 
minimum length of a set of vertex-disjoint paths that connect: 


x to ti , Sj to y, and Sh to th for every h = i + 1 ,..., j — 1 

These paths are illustrated in Figure [4} We interpret y as the last vertex of Pj[sj, bij] that is used 
in this sub-solution and we interpret x as the first vertex of Pj [bjj. tf\ that can be used in this 
sub-solution (or, more intuitively, the last vertex of the reverse of P t [b l _ 3 , t j]). By Lemma 3.1 each 
of the paths defining D t j[x, y] are contained by their respective region. 



Figure 4: An illustration of the paths defined by Djj [x;, y]. 


Optimal solution Given D\ we can compute the value of the optimal solution. By Lemma 3.4 
Q i and Qk contain a shortest (possibly trivial) path from si to a vertex x on Pi, and from a vertex 
y on Pfc to tk, respectively. Let y be the last vertex of Pk[sk,bi,k] that Qk contains and let x be 
the first vertex of -Pi[6i fc,ti] that Q\ contains. Then, by Lemma 3.4, the optimal solution has 
length Di k[x,y\ + £(a(si, x)) + l(a(y,tk)) where a is a Type I path between the given vertices. 
The optimal solution can be computed in 0(n 2 ) time by enumerating over all choices of x and y. 
Computing all such Type I paths takes 0(n 2 ) since there are 0(n ) such paths to compute, and 


each path can be found using the linear-time shortest paths algorithm for planar graphs HKRS97 


Base case: Leaf edges of T Consider a non-root leaf edge of T, which, by Observation 5.3 


is 


Bij -|-i for some i. Then Dij + i[x,y\ is the length of minimum vertex disjoint x'-to-f, and Si + \-to-y 
paths in Ri U i?,, + i. By Lemma |4.2[ Dij+i[x, y] can be computed in 0(n) time for any x and y and 
so D,jj + 1 can be populated in 0 (n 2 ) time. 


5.3 Non-base case of the dynamic program 

Consider a border B, j and consider the edges of T that are children of Bij. These edges considered 
counter-clockwise around their common node of T correspond to borders Bj j 1 , Bj 2 j 2 ,..., Bj t j t 
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where i < i\ < j\ < ■ ■ ■ < it < jt < j . For simplicity of notation, we additionally let jo = i and 
it +1 = j. Then, by Observation 5.2 either j# = ii + 1 or is the last border on Pj e and Bi e+1 ,j t+1 

is the first border on FT for l = 0 ,..., t. 


An acyclic graph H to piece together sub-solutions To populate Dij we create a directed 
acyclic graph H with sources corresponding to vertices of P, [bij , tf\ and sinks corresponding to 
Pj[sj, bij\. A source-to-sink (u-to-v) path in H will correspond one-to-one with vertex disjoint 
paths from: 


u to ti, Sj to v , and to t^ for every h = i + 1 ,..., j — 1 

Here u and v do not correspond to the vertices x and y that index Dij ; to these vertex disjoint 
paths, we will need to append vertex disjoint x-to-u and v-to-y paths (which can be found using a 
minimum cost flow computation by Lemma 4.2). 

The arcs of H are of two types: (a) Type I arcs and (b) sub-problem arcs. Directed paths in 
H alternate between these two types of arcs. The Type I arcs correspond to Type I paths and 
the endpoints of the Type I arcs correspond to the endpoints of the Type I paths. Sub-problem 
arcs correspond to the sub-solutions from the dynamic programming table and the endpoints of 
the sub-problem arcs correspond to the indices of the dynamic programming table (and so are the 
endpoints of the incomplete paths represented by the table). Note that vertices of a border may 
appear as either the first or second index to the dynamic programming table; in H, two copies 
of the border vertices are included so the endpoints of the resulting sub-solution arcs are distinct. 
Formally: 


• Type I arcs go from vertices of Pj e to vertices of Pi t+1 for i = 0,..., t. Consider regions Rj e 
and B.i t+ , • There are two cases depending on whether or not Rj t = Ri t+1 - 


— If Rj c = Ri e+1 , then for every vertex x of Pj e [sj e , bi e j e ] and every vertex y of Pj e [bi e+1 j e+1 , tj t ], 
we define a Type I arc from x to y with length equal to the length of the x-to-y Type I 
path. 

- If R jc / Ri e+1 , then for every vertex x of Pj e [sj e , b it j t ] and every vertex y of P ie+1 [bi e+u j t+1 , t je ], 
we define a Type I arc from x to y with length equal to the sum of the lengths of the 
x-to -tj e and Si e+1 -to-y Type I paths. 


• Sub-problem arcs go from vertices of P lf to vertices of P u for t = For every t = 

1 ,,t and every vertex x of Pi e and vertex y of P ]t (that are not duplicates of each other), 
we define a sub-problem arc from x to y with length equal to Di e j e [x,y\. 


Shortest paths in H By construction of H and the definition of Di e j e , for a source u and sink 
v , we have: 

Observation 5.4. There is a u-to-v path in H with length L if and only if there are vertex disjoint 
paths of total length L from u to ti, Sj to v, and Sh to t} r for every h = i + 1,... ,j — 1. 

See Figure [5] for an illustration of the paths in G that correspond to a source-to-sink path in 
H. Let H[u, v] denote the shortest u-to-v path in H (for a source u and a sink v). We will need to 
compute H[u, v\ for every pair of sources and sinks. Since every vertex in G appears at most twice 
in H, the size of H is 0(n 2 ) and for a given sink and for all sources, the shortest source-to-sink 
paths can be found in time linear in the size of P[ using dynamic programming. Repeating for all 
sinks results in an 0(n 3 ) running time to compute H[u,v\ for every pair of sources and sinksj^] 

2 Computing the length of the Type I paths is dominated by 0(n 3 ), but can be improved to 0(n log n) time by 


10 








Figure 5: Mutually disjoint walks represented by a directed path in H for a set of incident borders 
(green). The blue arcs correspond to the walks represented in sub-problems and the solid red paths 
correspond to the Type I paths represented by Type I arcs. The dotted red paths represent vertex 
disjoint u-to-x and v-to-y paths that will be added via a min-cost flow computation. 


Handling vertices that appear in more than two regions As indicated, a vertex c may 
appear in more than two regions; this occurs when two or more borders share an endpoint. In the 
construction above, if c appears in only two regions, then, c can only be used as the endpoint of 
two sub-paths (whose endpoints meet to form a part of an Sj-to-ij path in the global solution). 
However, suppose for example that c appears as an endpoint of both Bij and y and so 4 copies 
of c are included in H (two copies for each of these borders). On the other hand, one need only 
guess which Si-to-tj path c should belong to first and construct H accordingly. There are only k 
possibilities to try. 

Unfortunately, there may be 0(k) shared vertices among the borders 
Bi lt j 2 , Bi 2 j 2 ,..., Bi t j t involved in populating D iy It seems that for each of these 0(k) shared 
vertices, one would need to guess which Sj-to-fj path it belongs to, resulting in an exponential 
dependence on k. 

Here we recall the structure of T : the nodes of T correspond to non-regions', disks (or points) 
surrounded by regions. If there are several shared vertices among the borders, then there is an 
order of these vertices around the boundary of the non-region. That is, for a vertex shared by a 
set of borders, these borders must be contiguous subsets of B {, J2 , Bj 2J2 ,..., B lt j t . In terms of the 
construction of H , there is a contiguous set of levels that a given shared vertex appears in and 
distinct shared vertices participate in non-overlapping sets of levels. For one set of these levels, we 
can create different copies of the corresponding section of H. In each copy we modify the directed 
graph to reflect which s*-to-U path the corresponding shared vertex may belong to (see Figure [6]). 
As we have argued, since distinct shared vertices participate in non-overlapping sets of levels, this 
may safely be repeated for every shared vertex. The resulting graph has size 0(kn 2 ) since there are 
0(k ) borders and shared vertices are shared by borders. The resulting running time for computing 
all source-to-sink shortest paths in the resulting graph is then 0(fcn 3 ). 


Computing Dij from H To compute Dij[x,y\, we consider all possible u on Pi[x,U\ and v on 
Pj[sj,y\ and compute the minimum-length vertex disjoint u-to-x path and v-to-y path that only 
use vertices that are interior to Ri U Rj (that is vertices of Bij may be used); by Lemma 
paths can be computed in linear time. Let M[u,v\ be the cost of these paths. Then 


4.2 


these 


Dij [x, y\ = min M[u , u] + H[u, u]. 

uePi[x,ti],vePj[sj,y\ 


running Klein’s boundary shortest path algorithm 


Kle05 


in all regions, resulting in an 0(n 2 ) time to construct H. 
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Figure 6: (left) H constructed without handling the fact that vertices c\ and C 2 (black) may appear 
more than twice. The arcs are all directed upwards (arrows are not shown); green arcs are sub¬ 
problem arcs and red arcs are Type I arcs, (right) The levels that c\ appears in are duplicated, and 
only one pair of copies of c is kept in each copy of these levels. The vertex ci may only be visited 
twice now on a source to sink path. 


As there are 0(n 2 ) choices for u and v and M[u,v] can be computed in linear time, Dij[x,y\ can 
be computed in 0(n 3 ) time given that distances in H have been computed. 

Overall running time For each border Bij, H is constructed and shortest source-to-sink paths 
are computed in 0{kn 3 ) time. For each x,y G Dij[x,y] is computed in 0{n 3 ) time. Since 
there are 0(n 2 ) pairs of vertices in Bij, Dij is computed in 0(n 5 ) time (dominating the time to 
construct and compute shortest paths in H). Since there are O(k) borders (Lemma |3.3[ ), the overall 
time for the dynamic program is 0(kn 5 ). 

Acknowledgements This material is based upon work supported by the National Science Foun¬ 
dation under Grant No. CCF-1252833. 
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